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Non-topological solitons (NTSs) are extended objects that arise in theories with un- 
broken global symmetries [1,2]. The simplest example of a NTS is the so-called Q-ball 
that can appear in a C/(l)-invariant theory with a single complex scalar field 4> that has 
non-linear self-interactions [1,3,4]. The results for Q-balls found in the literature may be 
applicable to models that possess an unbroken global 1/(1) symmetry such as B — L. In 
this paper we study the properties of Q-balls in a local or gauged 17(1) theory. Our work is 
a first step in understanding how NTSs might arise in gauge theories such as electromag- 
netism, the Weinberg-Salam model, or Grand Unified Theories (e.g., SC/(5), SO(IO), E&). 

A NTS is a non-dissipative solution to the classical field equations that, for fixed 
charge Q, represents the field configuration with the lowest energy. For the Q-ball, the 
basic structure is very simple. Outside the Q-ball, C/(l) is unbroken and (f> particles have 
mass fi (this is the true vacuum of the theory). Inside the Q-ball, U( 1) is broken and a 
condensate of the 4> field forms with an energy per unit charge less than fi. Furthermore, 
the time-dependence of the 4> field gives rise to a non-zero charge density. The negative 
pressure of the false vacuum interior is balanced by the positive pressure from the confined 
<f> particles (for possible creation mechanisms see [5]). 

Coupling of <f> to a U( 1) gauge field leads to two basic effects. 1) Inside the Q-ball, the 
t/(l) symmetry is broken and the gauge field has a mass my. For Q-balls that are large 
compared to the Compton wavelength of the gauge field (£ = my), the charge is pushed 
to the surface of the Q-ball and the characteristic width of the shell in which most of the 
charge resides is 0{rn.y l ). It is in this limit that the Q-ball becomes superconducting. For 
small Q-balls, the charge is only slightly enhanced near the surface. These results Me a 
clear signal that the Q-ball is behaving like a lump of superconducting matter (which it 
must since U( 1) is broken inside the Q-ball). 2) Gauged Q-balls can only be so big. 

For small enough values of the charge Q and gauge coupling e, the ‘electrostatic’ 
energy of the Q-ball [=0(e 2 Q J /i?)] is smaller than the other energies in the problem and 
therefore we expect stable gauged Q-balls to exist. Suppose we fix e but slowly increase Q 
by bringing chaxges (i.e., free <f> particles) to the surface of the Q-ball from very far away. 
For some value of Q(= Q mo *) the cost in Coulomb energy for adding an extra unit of 


1 


charge becomes greater than the energy gained by bringing it from the true vacuum to the 
false vacuum. From this point on it becomes energetically favorable to leave additional 
charges as free <f> particles. A corollary of this result is that for large enough values of e 
there are no stable Q-balls. As will be discussed below, very small Q-balls are unstable 
and evaporate into free particles. Therefore, stable Q-balls must have charge greater -than 
some value Q min . If Qmax < Qmin then there are no stable Q-balls. 

The fact that there is a maximum value for the allowed charge of a gauged Q-ball 
points to a fundamental difference with the ungauged (e = 0) case. In the e = 0 case, it 
is possible to have bulk Q-matter occupying an arbitrarily large region of space. This is 
impossible in the gauged case. 

An outline of the paper is as follows: In Section II we describe the particle physics 
setting for gauged Q-balls. We review the basic properties of ordinary (e = 0) Q-balls 
and derive some general properties of gauged Q-balls. In Section III, we study gauged 
Q-balls using trial functions for the fields. The calculations are carried out for thin 
shell Q-balls and ones in which the gauge coupling and charge are assumed to be small. 
[How small will be made explicit below.] Section IV is devoted to numerical results. 
The equations of motion are numerically integrated on computer and some solutions 
are obtained. When appropriate, these results are compared with the analytic results 
obtained in the preceeding sections. In Section V we consider solitons with fermions that 
are coupled to a 17(1) gauge held. Finally, in Section VI, we summarize and discuss our 
results. 


2 General Setting 


Consider a complex scalar field <f>(r,t) = f(r,t)exp(i6(r,t))/y/2 coupled to a 17(1) 
gauge field A M . The Lagrange density for the theory is 

£ = - eArf - U(f) - (2.1) 

where — d^A u — and U(f) is the f7(l)-invariant scalar potential. [For definite- 

ness, we take e > 0.] The conserved charge associated with the J7(l) symmetry is 

Q = J d 2 rf 2 {6 - eA 0 ) (2.2) 
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where the dot means d/dt. For definiteness, we shall take Q to be positive. 

We begin by reviewing the basic properties of ordinary (e = 0) Q-balls. We assume 
that in the true vacuum f/(l) is unbroken, so that the absolute minimum of the potential 
[/(/) occurs at / = 0. The mass fi of a free (f> particle is given by 



(2.3) 


A necessary condition [4] for the existence of Q-balls is that 


= min 



(2.4) 


for some / = F a with F 0 finite and non-zero and 0 < u > a < fi. For large enough Q, the 
lowest energy state is a coherent configuration of the boson field <f> confined to a finite 
spatial region. Following Coleman[4], we refer to this object as a Q-ball. 

In order for a Q-ball to be stable against dispersion into free particles, its energy must 
be less than the energy of Q free particles. Consider a field configuration where inside 
some volume V, / has a constant value F and 9 = ut + 0 O where 6 0 is a disposable 
constant. Outside of V, f = 0. For large enough V (or Q) we can neglect the surface 
energies and treat / as a step function at the boundary of V. The total conserved charge, 
which characterizes the system, is Q = u}F 2 V. The total energy is then 


E = 


Q 2 

2 F 2 V 


+ U{F)V. 


(2.5) 


We minimize the energy with respect to V and find that 

V = i = (2.6) 

sJ2FHJ{F) 

and so the energy becomes 

* = (-) 

Next, we minimize E with respect to F keeping Q fixed and find that F — F 0 (cf. 
Eqn.(2.4)). For F = F a , the energy of the Q-ball is E = u> a Q. With u> 0 < fi, the Q- 
ball has lower energy than the energy of Q free particles and is therefore stable against 
dispersion. The negative vacuum pressure —U(F 0 ) is balanced by the positive ‘kinetic’ 
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pressure, u 0 F*/2. Assuming that there is no further complexity in the potential (e.g 
coupling of <f> to light fermions [6]), Q-balls will be absolutely stable. 

The result E = w 0 Q is valid only for very large Q-balls. For small Q-balls where 
surface energies are important, E can be much larger than u ) C Q and can in fact be larger 
than /aQ. In general there is a lower bound on the allowed values of Q for stable Q-balls. 
However, even for a large Q-ball, the surface energy is important as it determines the 
shape of the configuration. Since the surface energy is minimized by choosing, for a fixed 
volume, the shape with the smallest surface area, the groundstate Q-ball is a sphere. 

We now consider gauged U( 1) Q-balls (e / 0) where <f> is coupled to the gauge field 

as in Eqn.(2.l). One expects gauged Q-balls to be stable as long as their electrostatic 
self-energy is much smaller than the other energies. As before, / is non-zero and U( 1) is 
broken inside the Q-ball. As a consequence of the symmetry breaking, the gauge field is 
massive inside and the Q-ball is a U( 1) superconductor. 

Consider a coherent configuration of </> and A M with a given electric charge eQ. The 
lowest energy state will have no ‘electric’ currents and therefore no magnetic field. Fur- 
thermore, the lowest energy configuration will be spherically symmetric and stationary. 
For this configuration one can choose a gauge such that 9 = u?t ( 0 o — 0) with u ; constant 
and Ao(t') — ► 0 as r — » oo where r = |f|. For definiteness, we assume that u> > 0. The 
spatial components of the gauge potential are zero as there is no magnetic field. With 
the exception of 0, the fields are time independent. 

The Lagrangian for the configuration described above is 

L = in j rVr + i/V - 0(f)] (2.8) 

where y(r) = u ; — eA 0 (r) and prime denotes d/dr . By varying L with respect to / and g 
at fixed u? we find the equations of motion: 

f" + -/' + fg 1 - dU(f)/df = 0 (2.9) 

T 

g”+-g’-e 2 fg = 0. (2.10) 

T 

The total energy and charge are given by 

E = 4tt y r 2 dr ^ / ' 2 + ^ g 12 + ^ f 2 g 2 + U(f ) (2.11) 
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Q = 47r J r 2 drf 2 g (2-12) 

Eqn.(2.10) was derived by varying L with respect to Ao- Since Ao does not appear in the 
Lagrangian, we must treat Eqn.(2.10) as a constraint equation rather than a dynamical 
one. Eqn.(2.9) can be obtained by varying the energy in Eqn.(2.11) with respect to / 
under the constraint Eqn.(2.10). 

Some qualitative features of the soliton solutions are readily apparent. From Eqns.(2.10, 
2.12) we see that 

e 2 Q = lim 47 t r 2 g'. (2.13) 

For large r, g —+ u — e 2 Q/47rr, f is small and U(f) ~ /x 2 / 2 / 2. Eqn.(2.9) then becomes 

/" + -!' + (»’ - M 2 )/ = 0 (2.14) 

T 

and / oc exp( — T\J \i 2 — u? 2 )/r. For the existence of a localized solution without oscillations 
we require u ) < ji. Additionally, for the solution to be well behaved at the origin, /'(r) 
and g\r ) should approach zero at least faster than r for r — ► 0. 

We now show that g(r) obeys the following inequalities: 


0 < 5(0) < g(r) < g(o o) = u? < fi. 


(2.15) 


It is convenient to write Eqn.(2.10) in the following form: 



(2.16) 


We will first prove that 5(0) > 0. Suppose that #(0) = 0, Since y'(0) = 0, then g(r) 
would be zero for all r. Suppose instead that ^(0) < 0. Eqn.(2.10) then implies that r 2 g f 
is a decreasing function of r so that g'{r ) goes negative and g{r) < 0 for ail r. Neither of 
these possibilities are acceptable given that u > 0 and g(r) — * u) for r — ► 00. The only 
acceptable possibility is that 5(0) > 0. We then see that g f {r) is positive and therefore 
g{r) is a monotonically increasing function of r. 

Let us consider a different form of the energy integral. After partial differentiation 
and using the asymptotic behavior of / and 5, Eqns.(2.11,2.12) lead to the new energy 
functional 



(2.17) 
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where the constraint Eqn.(2.10) is satisfied. Eqn.(2.9) need not be satisfied. However, 
the energy is minimized by choosing / to be a solution of Eqn.(2.9). 

Demanding that the Lagrangian (2.8) for the solution is stationary at k = 1 under 
scaling of the form r — » kt leads to the following relation 

I dTT' [-i r + + 3 / [i/V - (/(/)] = 0. (2.18) 

When used together with Eqns.(2.9-2.12) and the asymptotic forms for / and g , the 
energy can be rewritten as 

E = wQ + y j dvr 2 [f 2 - g n ] . (2.19) 

Consider the case e = 0, so that g 1 — 0* In the thin wall approximation we can neglect 
the gradient f' 2 term. We then see that E = ujQ . It follows that for thin shell Q-balls 
with e / 0 the energy is generally bounded from above by wQ and is less than fiQ if 

U) < fl. 

Recall that our goal is to find the lowest energy configuration for a given value of the 
charge Q. One possibility is to have Q free particles of mass /a, the total configuration 
having mass fiQ. A second possibility is to have a soliton with charge 0 . In this case, the 
energy is a very complicated function of the charge. For e = 0 we know that E — ► u; 0 Q 
for large Q so that for < /x this configuration is favored over the free particle one. 
For e 0, we expect that the energy will be increased over the e = 0 case due to 
Coulomb repulsion with Coulomb energy becoming more important as Q gets large. For 
dE/dQ > fi we must consider a third possibility, namely that some of the charge can be 
put into the Q-ball and some can be put in free particles. Suppose that there exists a Q 
(= Qmax) such that for Q > £? mox , dE/dQ > fi while dE/dQ < fi for Q < Qmax • For 
Q > Qmax > a Q-ball with charge Q max plus Q — Q ma x free particles will be the lowest 
energy state for the system. 

As discussed above, in very small Q-balls surface energies are important and the Q- 
balls are generally unstable and disperse into free particles. There is therefore a minimum 
value allowed for the charge which we define to be Qmin • If the gauge coupling is large 
enough, Qmax will be less than Q m m> and there will be no stable Q-balls 
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Clearly, in order to understand gauged Q-balls, we need to determine E as a function of 
Q. We can do this analytically, but only in the thin wall approximation (i.e., neglecting 
surface energies) and for small values of e and Q. On the other hand, we can obtain 
better results numerically. The results obtained from these two approaches are presented 
in Sections III and IV respectively. 


3 Thin Wall Approximation 


In this section, we study properties of gauged Q-balls by considering tried functions for 
the fields. More specifically, we choose a simple trial function for the field / and solve for 
Aq (or equivalently, g ) using the constraint Eqn.(2.10). These functions can be thought 
of as an initial configuration for the Q-ball, and will in general have greater energy than 
the true groundstate configuration. [Recall that for the true groundstate configuration, 
the fields satisfy their equations of motion, Eqns.(2.9,2.10).] As we shall see in the next 
section, for a certain range of parameters, the trial functions used here closely approximate 
the actual solutions. 

We choose / to have a constant value F inside a sphere of radius R and zero outside. 
In real solitons, / changes continuously from F to 0 within a shell of some finite width 
T. For large Q-balls (width of the shell much less than the radius of the Q-ball) the 
energy associated with the shell is negligible (the shell energy is of order T/R of the 
total energy). In this limit we can treat / as a step function and discard the f n term in 
the energy density. This is the thin wall approximation. Aq or equivalently g must be 
determined by solving the constraint Eqn.(2.10). We find that 

/ \ [w — e 2 Q/47rR] R sinh(eFY)/r sinh(eFR) r<R , . 

^ r \ [a; — e 2 Q/4i:r] r > R 

where the gauge is chosen so that A 0 — ► 0 for r — ♦ oo. The system is now determined by 
four parameters, u Q,F, and R . In the following, we will determine the values for a ;,F, 
and R that minimize the energy given a fixed Q . 

Eqn.(2.12) for the charge gives one relation among these four parameters: 


u = 


e 2 Q 

47riZ 


1 - 


tanh x 


-l 


(3.2) 
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where x = eFR. We can use this expression to eliminate from E in Eqn.(2.17). 

i -i 


E = 


e 2 Q 2 
8tt R 


1 


tanh x 


+ y R 2 u(F) 


(3.3) 


Minimizing the energy with respect to R (or x) for fixed Q and F gives an expression for 
x in terms of Q and F: 

r _ -i 172 

(3.4) 


1 l] 

e 3 QF 2 

.tanhx 

4*/2 U(F) 


In principle, one can use this expression to eliminate R (and x) from the expression for 
the energy. The result would be an expression for the energy in terms of Q and F. The 
final steps are to minimize the energy with respect to F for fixed Q and to then eliminate 
F. The result is the desired expression for E in terms of Q. 

Unfortunately, Eqn.(3.4) is a transcendental equation for R in terms of Q and F and 
exact results from this point on are difficult if not impossible to obtain. We can however, 
carry out the above steps for the case where i < 1. As we shall see, x <C 1 implies 
that the Coulomb energy is smadl compared to the potential and kinetic energies. In this 
regime, the radius and total energy of the groundstate Q-ball differ from the e = 0 case 
by terms of order x 2 . We now proceed to calculate these corrections. 

For small x, Eqn.(3.4) becomes 

,2' 


1 - — 1 = e 3 QC 
15/ * 


or 


R = 


3 Q ) 


1/3 


^Fy/2U{F)) 


i + e 2 Q 2 ' 3 C 2 ' 3 ^ 


45 


(3.5) 


(3.6) 


where C = 3E'(E’ 2 /2U(F)) 1 / J /47r. For most of the potentials that we will consider, C is 
of order unity. Our expansion is therefore valid for x 2 . Note that when e = 0 we recover 
the result for the radius of ordinary Q-balls (Eqn.(2.6) with V = AvR 3 /3). 

We can use Eqns.(3.5,3.6) to eliminate x and R from Eqn.(3.3) to obtain 


E = Q i 


!2U(F) 


1 + 


e 2 02 /3 C 2 /3\ 


I' 


(3.7) 


F 1 V" ' 5 

Again, when e is set to zero we recover the result for the energy of ordinary Q-balls, 
Eqn.(2.7). 
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Let us assume that the function U / is minimized by choosing F — F a (cf. Eqn.(2.4) 
in the case where e = 0). The energy in Eqn.(3.7) is minimized for F — E 0 (l — A) where 

-f(sr([A(¥)L.r 

In deriving Eqn.(3.8) we have made use of the fact that d(U/ F 7 )/dF = 0 when evaluated 
at F = F a . By substituting F = FV(1 — A) into Eqn.(3.7) we obtain the desired expression 
for E in terms of Q. Likewise, we can use Eqn.(3.8) together with Eqn.(3.6) to find 
R = R(Q)- Eqn.(3.7) can be put into the following suggestive form: 


B ^ 


(3.9) 

F 2 20tt R k j 

The first term is the energy of the Q-ball neglecting gauge interactions. The second term 
is the electrostatic self-energy of the C^-ball, with e 2 Q 2 /87?R being the electrostatic energy 
in the region r > R and e 2 Q 2 /AOnR the electrostatic energy for r < R. Also in the small 
x approximation the first term is the leading order term. The second term is 0(x 2 ). 

In order to better understand the above results, we consider the following form for the 
potential: 

\ 2/6 /4 ,, 2/2 

(3.10) 


A 2 / 8 P n 2 P 

U(f) = J — + 

{J1 6 fi 2 4 2 


where A is a dimensionless constant and fi is again the mass of a free <f> particle. We require 
that A 2 > 3/16 so that U(f) > 0 for all f 7^ 0. As discussed above, a necessary condition 
for the existence of Q-balls is that the function U(f)/ f 2 has a minimum for some non-zero 
/. [It is for this reason that we must consider sixth (or higher) order potentials.] This is 
always true for the above potential and the minimum occurs at F a = \/3 m/ 2A. Note that 
the existence of Q-balls does not depend on the metastable false vacuum which exists 
when A 2 < 1/4. 

For ordinary (e = 0) Q-balls with an effective potential for 4> described by Eqn.(3.10), 
we get 

E = n*Q (3.11) 

where a = (1 — 3/16A 2 ) ly/2 = w 0 //x. Indeed, this energy is less than fiQ, the energy of Q 
free <f> particles. The energy gained by bringing a single (f> particle from the true vacuum 
to the false vacuum is just fi [1 — a]. 
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Given the potential, we can readily evaluate the quantities discussed above. We find 


that 


A = 


2 / e 3 QX 2 a 2 


2/3 


15 \ 7 r 

The energy in terms of Q is given by 


E = fiaQ 

Finally, the radius is given by 


30 \ 2 / 3 


20 \ir\a) 


(3.12) 


(3.13) 


R = 


1/3 Q 1/3 


fia 


1/3 


l + *-(»** 

45 V Tra 


2/3' 


(3.14) 


For a given charge Q , we see that both the radius and the energy are larger in the 
gauged case than in the e = 0 case. 

Let us now determine Qmax ? the maximum allowable charge for stable gauged Q-balls. 
From Eqn.(3.13), we find that 

dE 


dQ 


— fia 


- mr 


(3.15) 


We solve for Q ma x by setting this equal to /x: 

Q max ~ e 3 a l/2 ^ “ a ^ /2 • 

Using Eqn.(3.14), we can rewrite this expression in a more transparent form: 

e 2 Qr 


47rR 


= A [1 - a] 


(3.16) 


(3.17) 


The left-hand side is just the cost in energy due to Coulomb-type forces for bringing 
a particle to the surface of the Q-ball from very far away. The right hand side is the 
energy gained (when the Coulomb energy is neglected) in bringing a particle from the 
true vacuum into the false vacuum Q-ball interior. As stated above, and verified by 
Eqn.(3.17), Q max defines the value of Q where it becomes energetically favorable to keep 
additional charges outside and far away from the Q-ball. 

Of course Eqns.(3.12-3.17) are valid only when x is small. Recall that to leading order 


3 3 e^QF ( .F 2 y /2 

X 4tt \2 U(F)J ' 


(3.18) 
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For the potential Eqn.(3.10), we have that 


i 3 = 


3y/3 e 3 Q 
87rAa 


(3.19) 


At Q = Qmax i * = [3(1 — a)/a] 1/ ^ 3 . It therefore follows that x < 1 for a > 3/4 or 
A 2 > 3/7. 

From Eqns. (3.16) and (3.17) we can get the maximum radius for this configuration, 


ft™. = — t/— . (3.20) 

fie \ a 

In order to interpret the breaking of the U( 1) symmetry inside the Q-ball as super- 
conductivity, we need to compare the Umax with the penetration length £ associated with 
the mass of the photon (or equivalently, the Compton wavelength of the gauge field in a 
region where the t/(l) symmetry is broken). We find that 


1 1 _ 2A 

my eF 0 y/3fie 

(3.21) 

By comparing this with Rmaxi we get, 


3 

^ T ^ Rmax • 

4 

(3.22) 


We see that the results in this section are valid only in the regime where the penetration 
length is greater than the radius. In this regime, the charge distribution is roughly uniform 
being only slightly enhanced near the surface of the Q-ball. 

We conclude this section by considering the potential Eqn.(3.10) with A 2 = 3/16. As 
we now show, in this case £ << Rmax and the Q-ball is superconducting. For this case, 
a = 0 and the “false” vacuum (F = F a = 2/x) is degenerate with the “true” vacuum 
(F = 0). The potential energy is therefore zero inside as well as outside of the Q-ball and 
is non-zero only in the shell in which F changes from 2 fi to 0. Let T be the thickness 
of this shell and R be the radius of the Q-ball. From Eqn.(2.1l) we find the following 
approximate expression for the energy: 


E ~ 


8irfi 2 R 2 e 2 Q 


+ 


87ri2 


1 + 


v 


where 


/ tanhx V/ 1 1 \ 

' \x — tanhx ) \xtanhx sinh 2 xJ 


(3.23) 


(3.24) 
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Eqn.(3.23) deserves some explanation. First, we note that each of the terms in this 
expression is accurate only up to numerical factors of order unity. The first term is the 
scalar gradient term where we have set f f = F 0 /L — 2 ii/L* The second and third terms 
are the Coulomb and f 2 g 2 terms respectively. The last term is the potential energy in 
the shell where the average energy density in the shell is taken to be /x 4 /4. 

Minimizing the energy with respect to T we find that T = 2\/2a S ubstituting this 
into Eqn.(3.23) we find an expression for E in terms of R and Q. Let us now assume that 
x > 1. We will check this assumption at the end of the calculation. For x » 1 we find 
that A(z) ~ x“ 3 and 

-zrw / i \ 

(3.25) 

(3.26) 

(3.27) 
by solving 

(3.28) 


£ - 4 vw*- + ££ (i + s) 


Minimizing this expression with respect to R we find that 

2 -'I' 

\L V 87r7 

and 

1 (e<?) 4/3 

H 2 n /6^1/3 

where we have only used the leading order terms. As before, we find Q r 
dEjdQ = ji . We obtain 

Qmax = 

For weak couplings this can be quite large. 

Eqns.(3.23-3.28) are valid only for x 1 and we must check that this is indeed the 
case. For Q = Q max , x = which is indeed greater than one for weak (e < 2“ 3 / 2 ) 

coupling. For stronger couplings or for smaller Q the large x approximation breaks down. 

Finally, we compare the size of these Q-balls to the penetration length. From Eqn.(3.26) 
we see that 

jR l Q=Qm “ = 4v^//e 2 ^ 3 ’ 29 ^ 

This is greater than the penetration length ( = 1 /(e/r) for weak (e < l/(4\/2)) coupling. 
It is in this limit that the interior of the Q-balls becomes superconducting. From Eqn.(3.1) 
it is easy to see that the charge density at the center of the Q-ball is (D(e -X ) of the charge 
density near the surface and the width T of the shell in which must of the charge resides 
is of order (eF) _1 ~ (mv)' 1 as expected for a lump of superconductor. 
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4 Numerical Results 

In this section we present numerical solutions of gauged Q-balls. These results are 
used to 1 ) verify the general properties of gauged Q-balls derived in Section II; 2 ) check 
the accuracy of the thin wall approximation; and 3) explore properties of gauged Q-balls 
that cannot be studied in the thin wall approximation. 

The thin wall approximation discussed above is applicable to a rather limited set of 
Q-balls. Specifically, the results of Section III are valid for Q-balls in which the surface 
energy is small compared to the volume energy (i.e., the radius of the Q-ball is large 
compared to the thickness of the shell within which / changes from F to 0 ). However, 
for large values of e Q-balls have a maximum size that is small, so surface energies are 
always important. Moreover, many of the results in Section III assume that x<l where 
x is given in Eqns. (3. 18,3.19). This, for example, implies that the equations for Qmax are 
valid only for certain values of A. Because of the limitations of the thin shell analysis 
(the only analytical approach that we found tractable) the numerical analysis for gauged 
Q-balls is particularly important. 

Consider the Lagrangian Eqn.( 2 . 8 ) with the potential given in Eqn.(3.10). We scale 
the dimensionful quantities /, 5 , and r by appropriate factors of fi. In order to keep the 
notation simple, we use the same symbols for the dimensionless quantities in this section 
as we did for the dimensionful quantities in the previous sections (e.g., r here is equal to 
r//z of Sections II and III). For physical quantities such as E, we write the factors of fi 
explicitly. 

From the Lagrangian and the potential we get the equations of motion 

;|;W) = -sV + (a ! / s -/ s + /) (4-x) 

(4 - 2 > 

with the charge and the energy given by 

Q = 47T J drr 2 gf 2 (4-3) 


E = 



jor +£<✓>■+!/ v + 




(4.4) 
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As discussed above, a finite energy Q-ball solution must satisfy f'( 0) = (/(O) = 0 and 
5(00) < 1. Furthermore, a Q-ball solution with minimum energy for a given charge will 
be non-oscillatory (zero node solution). 

To find acceptable Q-ball solutions we numerically integrate Eqns.( 4 . 1 , 4 . 2 ) from r = 0 
to r = 00 (or rather, to a point well outside the Q-ball, i.e., where / ~ 0 ). [Recall-that 
we required A 2 > 3/16 so that !/(/) > 0 for all r.] Eqns.( 4 . 1 , 4 . 2 ) are equivalent to a 
system of four first order differential equations. In order to find Q-ball solutions, we must 
determine the boundary conditions (i.e., the values of f,g and their first derivatives at 
r = 0). We know that /'( 0) = <7'(0) = 0. Fixing <?( 0 ) determines the charge. /( 0 ) is then 
adjusted to give a non-oscillatory solution with /(r) — ♦ 0 for r — ♦ 00. [For definiteness, 
we choose /( 0) > 0.] 

Let us be more specific. For a given choice of /( 0), /(r) will (1) grow without bound, 
( 2 ) become negative for some r, ( 3 ) change from a decreasing to an increasing function of 
r, or (4) go smoothly to 0 a r — * 00. Note that the gauged case differs from the ungauged 
case in that f{r ) can begin as an increasing function of r and then turn over and go 
smoothly to zero at infinity. 

In practice, a value of /( 0) is selected, and the solution is integrated out until it 
becomes clear which of the four types it belongs to. 1 or 2 type solutions indicate that 
/( 0) was too big while type 3 solutions indicate that /( 0) was too small. The value of 
/( 0) is appropriately adjusted and the process is repeated until the desired accuracy is 
obtained. The value of <7(00 ) must then be checked to make sure it is less than 1. In Fig.l 
we plot f(r) and y(r) as a function of r for a gauged Q-ball with e 2 = 0.01 and Q ~ 10 4 

For comparison, we show on the same plot f(r) and for an ordinary (e = 0) Q-ball of 
roughly the same charge. Note that p(r) is a monotonically increasing function of r as is 
expected from Eqn.( 2 . 15 ). In Fig. 2, we plot the charge density p(r) — 2 gf 2 as a function 
of r for the gauged and ordinary Q-balls in Fig.l. While charge is distributed uniformly 
inside the ordinary Q-ball, the charge is pushed towards the surface of the gauged Q-ball. 
Note also that the radius of the configuration has grown compare to the e = 0 case. 
In Fig. 3, we plot E/pQ as a function of Q for various values of e. The unstable regions 
discussed in the text, those where dE (dQ > p are now apparent. Note that as we increase 
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the value of e 2 , the instability occurs for smaller values of Q. We also see that the slope 
for the e = 0 case never changes sign, indicating the absence of an instability for ordinary 
Q-balls. In Fig.4 we plot f[r) and p(r) for a Q-ball configuration using A 2 = 0.64 and 
e = 0.01. The charge for the configuration is 88,921. For comparison, we show the f(r) 
and g(r) used in the thin shell approximation. This indicates that the analytic ansatz 
used is not to far off. Finally, as mentioned above, there exists a value of e(= Ccrit) above 
which Q max < Qmin For e > e„i t there are no stable Q-balls. In Fig. 5. we plot e„i t as a 
function of A. 


5 Gauged Nontopological Solitons with Fermions 


The work presented above can easily be extended in a number of ways. The one we 
chose was to add Fermions to the previous model. Here we study nontopological solitons 
that occur when charged fermions are trapped inside a region of false vacuum[7]. Consider 
a theory with a real scalar field <r and Dirac fermion field ip. The scalar potential is taken 
to be 



(5.1) 


where we assume that a = 0 is the true vacuum and a = <t„ is the false vacuum. The 
Lagrange density for ip is 


= iipj* 1 (5 m + ieAp) ip — m ^1 — ipip 


(5.2) 


so that ip has mass m in the true vacuum and is massless in the false vacuum. 

The case where e = 0 was considered by Lee and Pang[7] and we briefly review their 
results before turning to the gauged case. Suppose there are N ip particles trapped in a 
spherical region of false vacuum. The energy for the configuration is 

3tt ( 3 \ 2 / 3 IV 4 / 3 2tt 


E 


4 V2t r 


J 


R 


+ y/ur, 


(5.3) 


The first term is the fermi energy for the ip particles while the second term is the energy in 
the surface of the soliton. [For the potential Eqn.(5.1) the false vacuum is degenerate with 
the true vacuum and so there is no volume term in the expression for the total energy.] 
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For fixed N the energy is minimized by choosing 


R = 


9 / 3 \ 2/3 iV 4/3 1 


16 V2tt ) 


j 


1/3 


At this radius, the energy of the soliton is given by: 



(5.4) 


(5.5) 


As discussed in [7], since the exponent of N is less than 1, solitons with large enough N 
will have energy less than N free “0 particles and will therefore be stable against dispersion. 
More precisely, for N > (3 18 tt 5 /2 19 ) (/xcr 3 /m 3 ) 3 , E < Nm and the solitons would be stable 
against dispersion into free particles. 

Suppose now that the fermions are coupled to a t/(l) gauge field A M as in Eqn.(5.2). 
The energy will be then given by Eqn.(5.3) plus a Coulomb term /?e N / R where /? is 
a number of order 10 -1 that depends on the charge distribution of the fermions. [For 
a uniform charge distribution, (3 = 3/20ir.j Again, we keep N fixed and minimize the 
energy with respect to R. We find that 


R = 


3 1.44iy 4 / 3 + /3e 2 iV 2 
47 r furl 


1/3 


(5.6) 


(lA4N i/3 +Pe 2 N 2 ) V3 (5.7) 

For stable solitons to exist, E < mN for some range of values of N. Clearly, this depends 
on e, the ratio /i<r 3 /m 3 and (3. e and /icr 3 /m 3 are model dependent parameters. (3 on the 
other hand, depends on the structure of the soliton and must be determined by explicitly 
solving for the soliton configuration. However, the exact value of /? is not necessary for 
the present discussion. 

In Fig. 6 we show E/Nm against N for two representative choices of e and /i<r 3 /m 3 [one 
in which stable solitons exist and the other in which they do not] and take /3 = 0.048(~ 
3/207 t). We see that fermion NTSs can occur even when the fermions are coupled to a 
U(l) gauge field. As with gauged Q-balls there is a maximum charge. Furthermore, in 
certain theories (c.f. example 2 in fig. 2) stable NTSs do not exist at all 


and 
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6 Summary and Discussion 

Solitons, whether topological or nontopological, may be interesting for cosmology as 
well as particle physics. The cosmological implications of topological solitons such as do- 
main walls, cosmic strings and magnetic monopoles have been studied for some time. On 
the other hand, the cosmological significance of nontopological solitons has been consid- 
ered only recently. The discovery of a relic abundance of NTSs would be clear evidence 
for a phase transition in the early Universe and therefore have important implications for 
particle physics. Moreover, NTSs are potential candidates for the dark matter in galactic 
halos [5]. It is therefore of great interest to explore the possible particle physics theories 
in which NTSs might arise. 

NTSs occur in theories with a continuous symmetry and therefore a conserved Noether 
charge. Previous investigations of NTSs have, for the most part, concentrated on theories 
with global symmetries. However, many of the theories that we know of (or at least believe 
in) involve gauge or local symmetries. In this work, we have considered NTSs (Q-balls) 
in the simplest gauge theory; namely, one with a local U( 1) symmetry. 

For fixed charge, the NTS is the groundstate configuration. Stability against dispersion 
into free particles is demonstrated by computing the energy of a soliton solution of charge 
Q and comparing this with the energy of Q free particles. Stability against fission, evap- 
oration of charge from the surface [6,8], and gravitational collapse must also be checked 
[9]. The existence of stable NTSs therefore depends on particle physics phenomenology 
such as the shape of the scalar potential or the coupling of scalars to fermions. This is 
in contrast to topological solitons where existence and stability are determined from the 
symmetry breaking pattern of a given theory [e.g., stable vortices or cosmic strings occur 
in any theory in three spatial dimensions in which a U( 1) symmetry is spontaneously 
broken]. 

Coleman [4] (see also Rosen [3]) has demonstrated the existence of stable NTSs or 
Q-balls in a theory with a single complex scalar field and a global U(l) symmetry. We 
have considered a similar model but with a local U[ 1) symmetry and find that gauge 
interactions affect the structure and stability of the Q-ball in a physically intuitive way. 
For small gauged Q-balls ( R m^ 1 ) configuration is perturbed only slightly from 
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the e = 0 case with a small charge enhancement near the surface of the Q-ball. For larger 
Q-balls (R >> rriy 1 ) most of the charge is distributed at the surface. This is consistent 
with the fact that the interior of the Q-ball is superconducting. Furthermore, there is 
a maximum size for stable Q-balls which occurs once the Coulomb barrier becomes too 
large. Finally we recall an important difference with the e = 0 case; the fact that.it is 
impossible to have gauged Q-matter. 
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Figure Captions 

Fig.l Plot of the scalar field /( i*) and the gauge field <7( r ) = w — eAo as a function of r. 

The solid lines are for e 2 = 0.01 with Q = 11,119. The dashed lines are for e = 0 

with Q = 10,941. For both cases, we have set A 2 = 0.2. 

Fig.2 Plot of the charge density p(r ) = 2 gf 2 as a function of r for the two field configu- 
rations plotted in Fig.l. Notice that in the gauged case, the charge accumulates at 
the boundary of the Q-ball. 

Fig. 3 Plot of E/pQ as a function of Q with A 2 = 0.20 for various values of e. 

Fig. 4 Plot of f{r) and g(r) (solid lines) for a Q-ball configuration with A 2 = 0.64 and 

e 2 = 10 -4 . Q for the configuration shown is 88,921. The dashed lines are the 
corresponding /(r) and < 7 (r) used in the thin wall approximation. 

Fig.5 Plot of e^it as a function of A. As discussed in the text, for e > e^it there are no 
stable Q-balls as Q m ax is less than Q m in- 

Fig. 6 E/mN as a function of N for the fermion ball. 
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